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Abstract 

Boundary quantum field theory is investigated in the Lagrangian framework. Models 
are defined perturbatively around the Neumann boundary condition. The analyticity 
properties of the Green functions are analyzed: Landau equations, Cutkosky rules together 
with the Coleman-Norton interpretation are derived. Illustrative examples as well as 
argument for the equivalence with other perturbative expansions are presented. 

1 Introduction 



In this paper we study quantum field theory with boundary in perturbation theory. The 
boundary is a fiat hypersurface with a space-like normal vector, i.e. in appropriate inertial co- 
ordinates it is just given by constraining one coordinate to take a constant value (conveniently 
chosen to be zero). The motivation for such a study and most of the explicit examples come 
from 1 + 1 dimensional field theory, but we present the formalism for an arbitrary number of 
spacetime dimensions to show its generality. 

Our main motivation comes from consideration of integrable boundary QFT in 1 + 1 di- 
mensions, although integrability is not required for the general formalism to work. In such 
QFT, it is possible to construct exact scattering amplitudes (S matrices and reflection ampli- 
tudes) using the bootstrap procedure. The central idea of the bootstrap is that certain poles 
of scattering amplitudes give rise to new states (bulk particles or excited boundary states) 
in the theory, which must be treated as fundamental degrees of freedom in their own right 
("nuclear democracy"). 
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However, in the general case there are poles which are not assigned to any particle. In the 
case of bulk theories, bootstrap closure (or consistency) means that all poles are explained in 
terms of intermediate processes involving particles on-mass-shell, and that the residues can 
be calculated using the on-shell couplings between the particles. In case of poles associated 
to bound states, the intermediate state is a one-particle state; the more general processes are 
associated to so-called Coleman-Thun diagrams [1]. Consistency requires that there exists 
a consistent choice for the on-shell couplings such that the singular contributions at all the 
poles can be calculated by the Cutkosky rules (in practice one normally looks only at the most 
singular terms). 

For QFT with boundaries, the respective rules for Coleman-Thun diagrams were general- 
ized by analogy in [2, 3], but there has been no systematic derivation of these rules so far. In 
the bulk, the usual derivation of singularity positions (Landau equations) and Cutkosky rules 
uses the analytic properties of covariant perturbation theory [4, 5], although the Cutkosky 
rules can also be understood as a generalization of unitarity [4]. 

There have been numerous previous works using perturbation theory for boundary QFTs 
[6]-[14]. Unfortunately these expansions are usually rather difficult to use to explore the 
analytic structure, because the propagators do not have simple analytic properties if expressed 
in terms of the momentum (in some cases even that is impossible). Therefore in Section 2 we 
propose to use a perturbation expansion around the free (Neumann) boundary condition and 
a constant bulk field configuration. It is not at all evident that this approach is correct in the 
general case, and the main difficulty comes from the fact that the vacuum of the boundary QFT 
may be a non-constant field configuration. As we mentioned above, expanding around such 
a field configuration would make the investigation of the analytic structure hopeless. Using 
a toy model, we present evidence in appendix A that a re-summation of our perturbation 
expansion gives the correct results that are expected from the standard approaches and in 
particular the tadpole contributions reproduce the effect of the nontrivial background field. 
As a further support for the correctness of the perturbative expansion, in Appendix B we give 
a systematic derivation of one-loop counter-terms in boundary sine-Gordon theory that have 
been used in semiclassical calculations in the literature [15, 16, 17]. 

We also extend the perturbation theory by taking into consideration fields living on the 
boundary, which is instrumental for the derivation of the boundary Cutkosky rules. We define 
asymptotic states, and write down the generalization of LSZ reduction formulae for boundary 
QFT, which were previously derived in a restricted context by us [18]. 

In Section 3, using the perturbative formulation, we derive the boundary extension of Lan- 
dau equations and Coleman-Norton interpretation, while in Section 4 we obtain the boundary 
Cutkosky rules. In Section 5 we present some explicit examples of singularities in the conjec- 
tured scattering amplitudes of boundary sine-Gordon theory. The paper ends by a discussion 
of the results and the outline of some issues that remain to be investigated. 

2 Perturbation theory with a boundary 

In the course of analyzing of boundary quantum field theories we follow the same line as 
presented in [5] for bulk theories. Setting the stage by introducing our conventions, we identify 
and canonically quantize the free theory. The R matrix is defined via asymptotic states and 
is related to the Green functions using the boundary analogue of the LSZ formula. Then 
we introduce Feynman rules for the computation of Green functions, both in coordinate and 
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momentum space. 
2.1 Conventions 

The coordinates describing the half space-time are 

z = (t, x, y) , x = (x 1 , . . . , x D1 ) , — oo < t, x 1 < oo, — oo < y < 
We denote the boundary coordinates by x = (t, x) and use the following abbreviations 



D-l 



dx; 



/«=n/ 

j i= i j- 



V = (d x i,...,d x D-i) 



dx 



oo 

dx I dt 

oo 



dz 



dx I dy 



Bulk fields are denoted by <& a (z), while for boundary fields we use (/> a (x). For simplicity, all 
the fields are supposed to be real scalars (generalizing to other cases is straightforward). The 
action is 

S = jdz\^- (d t <S> a ) 2 - (W a ) 2 - {d y <S> a ) 2 - Ml& a -F($ a )| 



U{<p a ^ a (y = 0)) 



(2.1) 



where 



M>3 {a lv . .a M } 

describes the bulk interaction, while 



« m 



u($ a ,4> a ) = Y, E E i; 

M,N {ai,...a M } {ai,...a N } 



ai...aM;ai...a^f^ai ■ ■ ■ ^om^cli ■ ■ ■ 0aj\r 



contains the bulk-boundary and pure boundary interaction terms (we suppose that it contains 
no terms with M = and N = 1, 2). 



2.2 Free fields 

Here we give a short description of the free theory with U = V = 0. Free bulk fields satisfy 
the equations of motion 

(d 2 -V 2 -d 2 y +M 2 )<f> a (z) = , d y * a (z)\ y=o = 

The field satisfying the boundary condition can be decomposed as 

$ a (t,X,y) = J ^ COs(/q/) J ^D-l ex P(^ ' X)$ a {K, k, t) 
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Similar decomposition is valid for the canonical momentum H a (z) = dt&(z). The commuta- 
tion relations are 

[$ a (t, x, y),U p (t, x, y')} = i5 a(3 5(x - x") (5(y - y') + 5(y + y')) 

due to the Neumann boundary condition. We can introduce the bulk creation/annihilation 
operators 

A a (n,k,t) = iIL a (n, k,t) + Q a (n, k)& a (n, k,t) 
A a (K,k,t)^ = —ifl a (K,k,t) + Q a (K,k)^ a (K,k,t) 

that satisfy 

A a (K,k,t), Ap(K',k',tf] = (2tt) d 2tt a (K,k)5 afS 5(k - k') (5(k - k') + 6(k + k')) 

The boundary fields can be quantized in the usual way as free fields living in D dimensional 
spacetime: 



4>a{x) = 



dk 



(27T 



^-exp^/c ':■ x)(f) a (k,t) , 4> a (k,tf = 4> a (-k,t) 



= 9 t (j>a , [<f>a(t, x),TT b (t, x')] = iS ab S(x - x") 

a b (k,t) = m b {k,t) +Ub(k)4>b(k,t) , uj b (k) 



k 2 + m 2 



a b (k,t), o c (P,t)t] = (2n) D ' 1 2co b (k)S bc S{k - k') 
The (normal ordered) free Hamiltonian can be written as 

H= \j 2%\ (2tt) d ~ 1 Aa(yK ' ^ *) ti4a (*' ^ + / ^D-i ab $itf a b{k,t) 

and the time development of the modes can be calculated as 

A a (K, k, t) = e-^^A^K, k) , a b (k, t) = e'^^a^k) 

The Fock space of free fields can be introduced in the standard way. The vacuum satisfies 

A a ( K ,k)\0) = = a b (k)\0) 

and the space is spanned by the states 

A ai (m, fci) f . . . A aM (K M , k M y 'a bl (k[y . . . a bN (k' N ) ] |0) 

Due to the symmetry A a (K, k) = A a {—K, k) we could constrain k > 0. However, it turns out 
to be simpler to let k take general real values, and impose the former symmetry property to 
account for the presence of the boundary. 
The free propagators are 



9ab{x,x') = (0| Tcj) a (x)(f) b (x') |0) = i5 ab 



dk e lfc- ( x x ') 
(2ir) D k 2 -m 2 + is 



where k = (ko, k), J dk = J dko J dk, and 



G af} {z,z') = {0\T<f> a (z)<!>p(z')\0) 

- i5 „ r — l ^ e^^2 ( e i K (y-y>) iK (y + y>)) 

~ %daf, J_ 00 2irJ (2*)°k*-K*-Ml + ieV + ) 

The bulk propagator contains two terms: the first one describes the direct propagation of 
the fields, while the other corresponds to reflection on the free boundary. This can also be 
interpreted as a propagation to the mirror image point at — y', and will play an important role 
in the formulation of the Feynman rules later. 

2.3 Asymptotic states, R matrix and reduction formulae 

We proceed with the usual assumptions: the interaction is switched off adiabatically as t — > 
±00. The asymptotic fields are then defined to be free fields and 4>^ s (as = in, out), 
satisfying 

lim ^ a (t,x,y)-Z a ^ out (t,x,y) = 
lim M^^-Z^^^x) = 

where the Z coefficients, as usual, take care of the canonical normalization of the fields, and 
the limits are understood in the weak sense (i.e. for any matrix element of the fields). 

The modes of the asymptotic fields, which are obtained using the free field mode expansions 

A a a s ( K ,k) = 2i f° dy f dxcosinyy^^-^Vt^ (t,x,y) 

J —CO J 

A a a s (K,k)^ = -2% f° dy f dxcos^e-^^ 1 - 1 ^^ 3 (t,x,y) 

J —CO J 

a a b s (k) = i J dxe l ^ t ~ t£ ^ t <p a b s {t,x) 

create asymptotic states as follows: 

m, ki, ai; . . . ; k m , k M , olm; h, h; ...;k N , b N ) = 

I as 

AZ ( Kl , ^)t . . . A% s Jk m , k M )^all{k^ . . . a a b s N {k N y |0) . 

We assume asymptotic completeness: both the in and the out states form a complete basis. 
The unitary transformation between the two is what we call the reflection matrix R (in the 
usual terminology, reflection matrix means the matrix element of R between a bulk one-particle 
in and a bulk one-particle out state). For any given initial and final state, the corresponding 
matrix element of the R matrix (more precisely: R operator) gives the probability amplitude 
for the evolution of the initial state \i) at t = — 00 into the final state |/) at t = +00: 

Rfi = out {f I i)in 
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Now one can proceed to deduce reduction formulae for the matrix elements between asymptotic 
states. The derivation is essentially the same as in [18], so we only write down two examples 
of the resulting formulae. Applying the formalism to an incoming bulk particle gives 

out (A\0\ K ,k,a;B) in = out {A\OA%( K ,k?\B) in = 

= disconnected part + 2iZ~ 1/2 / dy dx cos{Ky)e~ iiSla{ ^ k)t ' k ^ ) x 

J —CO J 

[d 2 - d 2 y -V 2 + M 2 + 6(y)d y } out (A\ TO$ a (t,x,y)\B) in 

where O denotes a general T product of local operators. The reduction formula for an incoming 
boundary particle created by the a modes are identical to the usual reduction formulae in D 
space-time dimensions without boundary: 

out (A\0\k,b;B) m = out (A\Oai n (k)^\B) in = 

= disconnected part + iZ' 1 / 2 j d xe -i{^ik)t-k-3) x 

[d 2 - V 2 + m 2 a } out (A\TOMt,x)\B) m 

Formulae for outgoing particles can be written in a very similar form. 

2.4 Feynman rules in coordinate space 

In the interaction picture, the time evolution of the system can be described by the operator 
U(t) = Texp^-iJ dTH int (T)^ 
H mt (r) = £ dy J dxV (<&™(r, x, y)) + J dxU (<&™(r, x,y = 0), x)) 
The R matrix can be expressed as 

R = U(oo) = T exp | -i J d T H int (r) j 

and the interacting fields take the form 

$ a (t,x,y) = Uitr^Zfr^yMt) 

<t> a (t,x) = u(t)- i 4?:(t,x)u{t) 

Using the free field formulae and Wick's theorem, one can readily derive the Feynman rules 
in coordinate space for the Green's functions 

G a l...a m ;a 1 ...a n ( x l> Uli^l ■ ■ ■ i x mi Dmi ^mj x li tli ■ ■ ■ x ni ^n) = 

(0| T$ Ql (fi,yi,ti) . . . <f>a m {Xm,ym,t m )<l)a 1 (x'l,t' 1 ) ■ ■ ■ 4>a n {x' n ,t' n ) |0) 

from which the rules for the R matrix can be obtained by applying the reduction formulae. 
The resulting diagrams contain the following ingredients: 
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1. Boundary propagator: 

a b f dk e -ik-(x-x') 

t, x t',x' = l ab J (2tt) d k 2 - ml + is 

2. Direct bulk propagator 

a _ ? = „ f° - /" d * 

f,y,a? aP J. 0O 2irJ (2ir) D k 2 -K 2 -M* + ie 

3. Reflected bulk propagator 



/3 [~dK[ dk e-tt-C-*') ^ 

f,y,£ ^J-oo^J (2ir) D k 2 -K 2 -ml + ie 

4. Bulk vertex (it doesn't make a difference whether a leg comes from a direct or reflected 
bulk propagator) 

.C*3 




Qf4 f ai...«M 

5. Boundary vertex (once again, no distinction between the types of bulk legs) 

'. ' = ~ iu a 1 ...a M ;a 1 ...a N 
OLM Q>N 

The vertex positions must be integrated over, for bulk vertices over the bulk (half) space, 
while for boundary vertices over the boundary. There are also combinatorial factors that 
follow straightforwardly from counting the number of ways the particular Wick contractions 
can be made. One has to take care to draw all possible diagrams obtained by putting both 
the direct and the reflected propagators for each bulk leg or line. In fact, for coordinate space 
rules we could have used one type of bulk propagator as well: G a p(z,z'), which consists of 
the sum of the two pieces. The separation of the bulk propagator into two pieces, however, 
makes the momentum space formulation much easier. 

2.5 Feynman rules in momentum space 

We can give the Feynman rules in momentum space by taking a simple Fourier transform of 
the Green's functions: 

G™{ n .. am .-,ai...a n (^!' Kl > UJl ■ ■ ■ ' Km ' L ° m '^ k[,uJi, ■ ■ ■ k' n , u!^ = 

m / p+co P \ J 1 / P 

n 1 / d Vi dudx.s^-^-^ ) n ( / dtjdxjjWi-^ 

i=\ J— 1 

Gai...a m ;a 1 ...a n { x liVlitl ■ ■ ■ j ii Urn, tm]Xi,t 1 , . . . X n , t n ) 
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(the y integrals are denned by a simple extension of the Green's functions as even functions 
of y). 



1. Boundary propagator: 

a b _ i8 ab 

u>,k to 2 — k 2 — m\ + is 

2. Direct bulk propagator 

a P i8 aj3 



ft, K ,k tt 2 -k 2 -K 2 -M 2 + is 

3. Reflected bulk propagator 

a (3 i5 a /3 

ft, k, k " ft 2 - k 2 - k 2 - Ml + is 



4. Bulk vertex 

a 4 = -tv ai ... aM (2ir) D 5( D \Zk)ir5(E' «) 

The momentum conservation for the D-momenta is the usual one, but for the n 
component, the prime indicates that if two vertices are connected by a reflected bulk 
propagator, the corresponding k must be oriented either outgoing or incoming at both 
vertices. For direct bulk propagators, momentum is conserved as usual: all D + 1 
components of the momentum are oriented as outgoing at one end and incoming at the 
other one. 

5. Boundary vertex 

///,,..., l ,, : „,.,, A (^) / v / Nv/,-) 

OLm 0-N 

There is no conservation law for the component k, since the position of boundary vertices 
is only integrated over along the boundary. 

There is also an integration over the momentum of each internal line. 

As a result of the above rules, the presence of reflected bulk propagators and boundary 
vertices breaks k conservation, as expected from the absence of translational invariance in the 
direction perpendicular to the boundary. 



3 Landau equations 

The analyticity properties of the Green functions can be analyzed in the perturbative frame- 
work. The Landau equations describe the singularity structure of the various Feynman dia- 
grams which we are going to derive for the boundary theory. 
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3.1 Derivation of the Landau equations 



Consider a Feynman diagram in our generic scalar theory (2.1) with N outer bulk legs with 
momenta (ki,Ki), i = 1,...,N and M outer boundary legs with momenta hi , i = N + 
1, . . . ,M + N, where fcj are the components of the momenta parallel to the boundary and 
Ki are the transverse components (only for bulk legs). Using the fact that D-dimensional 
Poincare invariance (parallel to the boundary) is unbroken, the resulting amplitude depends 
only on the invariants ki ■ kj and Kj. The general Feynman integral can be written 

. L D K , I I+J 

° = j n n # n^ 2 - - m? + ^r 1 n & -r 1 m 



=1 r=l s=I+l 



where % denote D- dimensional loop momenta for L loops in total, Xj are the transverse "loop" 
momenta (see later), (p r , 7r r ) are the D + 1-momenta of I internal bulk lines, while p s are the 
D-momenta of J internal boundary lines. Momentum conservation can be used to expressed 
p and 7r in terms of q, x, k and k (taking into account that a reflected bulk propagator reverses 
transverse momentum). 

We can then introduce Feynman's parameterization in the usual way 

g = / n (Cp> n n jf (e - i ) «> «> * a )~ j ~ J 

J 7+ J 

4>{k,K,q,X,a>) = J2 a r(Pr -^r ~ M r + ie ) + E a s(Ps- m s+ ie ) 
r=l s=J+l 

and integrating out the loop momenta g, % (these integrals are - apart from UV divergences 
taken care by counter-terms - nonsingular in general, just as in the case of bulk diagrams), 
we are left with a multiple integral over the a, which can be considered as a function of the 
invariants ki • kj and Kj. 

Following the usual argument, a singularity can occur (in the limit e — > 0) when the hyper- 
contour in a space is trapped between two singularities of the integrand, or a singularity of 
the integrand takes place at the boundary of integration. The integrand is singular if 

a r = or p 2 — it 2 — M 2 = , r = 1, . . . , I 

a s = or p 2 s — M 2 = , s = I + 1, . . . , I + J (3.2) 

which are nothing other but straightforward generalizations of the usual Landau equations. 
Reduced diagrams are defined by shrinking internal lines with a% = to a point, and the 
remaining lines must be on the mass-shell. A further requirement is that 



d_ 

dqi 



/ I i+J \ 

\r=l s=I+l / 



and 



d ( 1 I+J ' 

— a r (p 2 r . - it 2 - M 2 + ie)+ a *(Ps ~ m l + ie ) 

° Xi \r=l 8=1+1 



Figure 3.1: Generic form of a closed loop 



These equations can be brought into a more elegant form. First concentrate on equation 
(3.3) that is on the Lorentz invariant part. Consider a closed loop in the Feynman diagram 
(Fig. 3.1) where the internal momenta are relabeled, the vertices can be of boundary or of 
bulk type and the various l-s collect all the outer momenta from the loop point of view, but in 
principle they can be outer (k) or inner (p) or their sums. Using momentum conservation the 
momentum integration can be evaluated. As a result every pi , i = 2 . . . A can be expressed in 
terms of p\ (which serves as the loop variable) as follows 

A-l 

P2 = pi + h ; P3 = P2 + fa = pi + fa + fa ; • • • pa = pi + ^ h ( 3 - 5 ) 

3=1 

and the overall momentum conservation Ylik = m ust hold. Substituting (3.5) into (3.3) for 
Qi = Pi we obtain 

^ a iPi = (3.6) 

each loop 

This argument is valid for a D dimensional bulk theory as well as for the D non-transverse 
component of the momenta in a boundary theory. 

Now lets focus on the transverse momentum, that is on (3.4). Consider a loop as before but 
with bulk vertices only (Fig. 3.2). For bulk propagators we take e = 1 while for the reflected 
one e = — 1. Using momentum conservation the momentum integration can be eliminated 
successively giving rise to 

7T2 = tlTTl + h ; VT3 = 627T2 + fa = ^2^1 + ^1 + fa ', 

and in general 

it i 

TTi = Y\ ejTTl + II tjh H h Yl e 3 l k ^ h l * 

3=1 3=2 j=k+l 

The delta function 5 (eA^A + I A — tti) gives tt± = ha^i + where we have introduced 
lii = ni=i € j to encode the parity of the momentum change caused by reflected propagators. 
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Figure 3.2: Generic form of a closed loop with bulk vertices only 




Figure 3.3: Path starting from and ending at a boundary vertex 

Clearly if ha = — 1, that is the loop contains odd number of reflected propagators the delta 
function eliminates the last integration determining 7Q in terms of the l-s. In this case there 
is no singularity. If, however, ha = 1 (loops with even number of reflected propagator) tti is 
not determined, we have to integrate over it (tti = xi) and we obtain the Landau equation: 

HiOiTTi = (3.7) 

each loop 

Now consider a path in the Feynman graph starting and ending with boundary vertices but 
containing bulk vertices only (Fig. 3.3). Using the previous calculation tta can be expressed 
in terms of 7Tias before tta = HA-i^i + • • • but now there is no more delta function so tt\ serves 
as a "loop" variable, %i, for which integration has to be performed. As a result we obtain the 
following Landau equation: 

^ HiOLiTii = (3.8) 

each path 

There is an interpretation of equation (3.6) in the bulk theory in terms of electric network: 
it is just translated to the usual Kirchoff laws for the .D-momenta of the reduced diagram 
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with "resistances" a r and "currents" p r so the "voltage drops" are a r p r . The generalization 
for the boundary case: equation (3.7) means that a r -s are the "resistances" and the 7i>-s are 
the "currents". Note that there is one current for every component of the D + 1 dimensional 
energy-momentum vector and so for every loop and vertex there are D+l independent Kirchoff 
laws to write down. 

The signs introduced by the reflected propagators can be traced as follows: every reflected 
propagator carries an "inverter" device that changes the sign of the transverse momentum 
current, so fi r keeps track of the actual sign of the current and fi r a r ir r are the "voltage 
drops". Equation (3.8) means that as far as the transverse momentum current is concerned, 
all the boundary vertices are on an equipotential surface, which can be considered as the 
reference point ("earth"). For the components of the energy-momentum that are parallel to 
the boundary, the Kirchoff laws take the usual form as in the bulk theory. 

3.2 Coleman-Norton interpretation 

The Landau equations for a generic Feynman graph in the bulk theory are those we have in 
equation (3.2) and (3.6). Their physical region singularities have «i > and contain real 
momenta pi. Following Coleman and Norton [19], these equations can be interpreted as the 
existence of a space-time graph of a process involving classical particles all on the mass shell, 
all moving forward in time, and interacting only through energy and momentum conserving 
interactions localized at space-time points. The correspondence is one-to-one and goes as 
follows: for each internal line draw a vector a r p r , of length a r M r . (Lines with a r = are 
shrunk to a point). For each vertex in the graph a space-time point is associated where the 
momentum conserving interaction occurs. The consistency of this picture means, that two 
different paths leading to the same vertex in the Feynman graph define the same space-time 
point, which is nothing but the equation (3.6). 

The Coleman-Norton type interpretation of the boundary Landau equations is the same 
as in the bulk case: the existence of a space-time graph of a process involving particles all 
on the mass shell, all moving forward in time, and interacting only through local interactions 
localized at space-time points where now particles can scatter on the boundary, which reverse 
the sign of the momenta. The correspondence goes as follows: For each internal bulk line draw 
a vector (a ) with length a r M r (reflected line in the case of a reflected propagator), 

while for boundary lines a vector (a s p s ) of length a s m s (lying in the boundary). Bulk vertices 
are located at space-time bulk points where the momentum preserving interaction occurs. 
Boundary vertices are located at the boundary points with interactions where the transverse 
component of the momentum is not conserved. To see that for a space-time diagram the 
Landau equations are satisfied consider the diagram in Fig. 3.4. 

On the left side of the boundary is the real space-time picture of some process, while on 
the right hand side is its mirror image with respect to the boundary, which is necessary to 
introduce in order to interpret it in terms of the Landau equations. Instead of drawing a 
reflected line for a reflected propagator (e.g. from 1 to 2) we consider a straight line from 
the real world to the mirror one (from 1 to 2') or vice versa. The variable // keeps track on 
which side of the boundary we are and to which "world" (real or mirror image) the next line 
is directed. Consistency of the Coleman-Norton picture means that for two different paths 
between two interaction point the displacement vectors must sum up to the same value, which 
is ensured by equation (3.7). Evidently we compare only points that are both on the "real" side 
of the boundary so it is sufficient to consider loops with even number of reflected propagators. 
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Figure 3.4: Space time diagram of a general on-shell process using the method of images 



Another consistency condition is that all the boundary vertices must have the same transverse 
coordinate, which translates into equation (3.8). 

In order to reverse the above argument, we must show that to any solution of the Landau 
equations, a space-time diagram can be associated. This can be done by drawing the extended 
diagram first, similarly as we would draw in a bulk theory, then choosing a loop with odd 
number of reflected propagator and determining the location of the boundary by identifying 
the endpoints of the closed loop as the mirror image of each other. Clearly such a loop is 
not closed in the transverse coordinate. Now the Landau equations for the loops with even 
number of reflected propagators ensure that the location of the boundary does not depend on 
the choice of the loop with odd number of reflected propagators. 

4 Cutkosky rules 
4.1 Derivation 

A further important ingredient is provided by the boundary analogue of Cutkosky rules to 
calculate discontinuities across the cuts in the physical region determined by the solutions 
of the Landau equations. To give the boundary generalization of Cutkosky rules we need a 
further definition. Note that if we write down a Feynman integral G of the form (3.1) with 
I bulk internal lines, each bulk line corresponds to either a direct or a reflected propagator 
so there are actually 2 1 distinct diagrams that differ only which bulk lines are taken to be 
direct or reflected (represented by continuous resp. dashed lines in the associated Feynman 
diagram). Let C{G) denote the sum of terms the possible assignments of direct /reflected 
bulk lines. Each individual term can be written in the form (3.1), the difference being the 
expression of the loop momenta \ m terms of the internal momenta k due to the different 
conservation laws obeyed by the two different types of bulk lines. 

It is clear that C(G) can have a singularity at a given value of external momenta only if 
some individual term in the sum does. For an integral of the form (3.1), every singularity 
corresponds to a solution of the Landau equations (3.2,3.3,3.4). In such a solution there are 
internal lines which have a corresponding a r that is nonzero (and so are on the mass shell); 
these are the ones that make up the reduced diagram corresponding to the given singularity. 
We remark that since the momentum conservation conditions (3.4) are different for each 
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individual term in the sum C(G), in general only a few terms give a nonzero contribution to 
a given singularity. 

Let us denote the set of lines which are on-mass-shell at the singularity by Go. We restrict 
ourselves to cases when the graph G\Go, obtained by severing the lines on-mass-shell separates 
into two disconnected graphs G\ and G2 (in the bulk the corresponding singularities are called 
"normal thresholds") . 

Then the rules are that the discontinuity of G(G) (3.1) across the cut corresponding to 
the singularity can be evaluated by substituting the propagators of the internal lines for which 
a r 7^ as follows: 

-2md^ (p 2 r - n 2 - Mf) 

p 2 r - 7lf _ M 2 + K \fr r r) 

-= \ -27ri5^ (p 2 r - m 2 r ) 

Pr — +ie 

change e — > — e for internal lines in one of the components (say G2, the other choice just gives 
the jump in the opposite direction across the cut), and then perform the integration in (3.1). 
Here 5^ denotes a delta function where only the root po = +\ff% + tt 2 . + M 2 or + \/p 2 . + m 2 
is taken into account. 

To show this, we generalize a very elegant proof for the bulk Cutkosky rules, originally 
due to Nakanishi. The bulk proof is spelled out in detail in [5], so we only give the details 
necessary for its extension to the boundary case. We construct an auxiliary field theory to 
any given diagram class C(G) in the following way. To each internal bulk line of the diagram 
we associate a different species of bulk field $ r of mass M r (r = 1, . . . , I), similarly to each 
internal boundary line a boundary field <j) s of mass m s (s = 1, . . . , J). To every bulk vertex 
v with some external bulk line we attach a field <& v of mass M v , where M 2 = P v is the 
total (D + l)-momentum entering the given vertex from its external legs. Similarly to every 
boundary vertex u we attach a boundary field <p u of appropriate mass m u , rfi^ = p\ where 
p u is the total D-momentum entering the vertex from external bulk and boundary lines. The 
auxiliary action corresponding to C(G) is defined as 

r=l v L v r^v 

I dx e \ [^) 2 - m2 ^] + e % \{ d ^) 2 - ^«^| + e ^ n *r n < 

s=l u u r— »u s— *u 



's 



where Y^' v an d Yl'u mean that the sum goes only for bulk/boundary vertices that have external 
lines entering them, and (p u must be put equal to 1 for vertices which connect only to 
internal lines, and expressions of the type r — > u means taking product over internal lines r 
that connect to vertex u. 

The crucial point is that C(G) gives the lowest order contribution to the scattering ampli- 
tude with the given external momenta in the auxiliary theory. Let i denote the initial state of 
this scattering and / the final one (external lines attached to Gi and G2, respectively). Then 
the transition amplitude for this process can be written 

T fi = C(G) 
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However, on-shell processes are unitary, thus we have 

Tf l -V s = Y,C(G 1 )C{G 2 y 

G 

where X^g means a phase space summation over the lines involved in Go- Following a similar 
reasoning as in [5] yields exactly the rules spelled out above. 

4.2 Comments on the general case 

For different topologies, we expect that similar rules hold. The reason is simple: the Cutkosky 
rules themselves depend only on the analytic structure of propagators. In the perturbation 
theory around the free (Neumann) boundary conditions we see that the momentum space 
form of propagators is just the same as in the bulk theory, the only difference comes from the 
momentum conservation rules which express the internal momenta in terms of the external 
and loop momenta and from the fact that boundary propagators have vanishing transverse 
momentum. However, insofar as these expressions are linear (which is true even in the bound- 
ary situation), neither the precise form of the conservation rules enters in any derivation of 
the Cutkosky rules, nor the constraint on the transverse momentum for boundary lines (which 
is also linear) is relevant. Therefore we can see that the boundary Cutkosky rules given above 
can be derived for any topology for which a derivation is given in the bulk (e.g. the well- 
known triangle anomalous threshold diagram, which is not in the class covered by Nakanishi's 
derivation) . 

In the above derivation we also supposed that and rh^ are positive, which corresponds 
to singularities in the physical region. Beyond the physical region nontrivial analytic continu- 
ation is required similarly to the bulk case. Extension of Cutkosky rules to such singularities 
is a very complicated issue, tantamount to a generalization of unitarity. In this paper we do 
not enter into further discussion of singularities in the non-physical region, but simply suppose 
that an extension of these rules can be worked out. 

Finally we note that singular contributions corresponding to the same reduced diagram 
can be summed up, giving the familiar result that the total contribution (to all orders of 
perturbation theory) can be obtained by substituting the exact vertex functions for the vertices 
in the reduced diagram characterizing the singularity (see e.g. Section 2.9. of [4]). 

One has to be very careful in deducing the singularities of the Green functions from the 
analyzis of the perturbative series. Sometimes individual terms may have singularities that 
are cancelled when a sum of terms is taken. Some other times individual terms do not yield 
the correct singularities at all. It happens when some nonperturbative dynamical property, 
such as formation of a bound state, is involved. To handle these situations a complete or a 
partial summation of the perturbative series has to be performed. 

4.3 Infrared divergences 

As shown in appendix A. 2, perturbation theory around the Neumann boundary condition has 
infrared divergences. These arise from on-shell bulk particles propagating along the boundary 
as expected from the Landau equations and correspond to a boundary bound state. 

However, if there is a quadratic term of the form \&(y = 0) 2 /2 for the corresponding 
bulk field $ in the boundary potential V (which is the case in general), then the infrared 
divergence is nonphysical. There are infinitely many infrared divergent terms with different 
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Figure 5.1: Diagram corresponding to an on-shell boundary state 

degree of singularity and resummation in the coupling A shifts the singularity from k = to 
k = —iX, and it does not affect the derivation of the physical region Cutkosky rules presented 
above. 

For A > the singularity is shifted to the unphysical sheet of the scattering amplitude, so 
it is physically irrelevant. For — m < A < 0, the pole is again shifted away from the region 
of physical values of the incoming momentum, however, it remains on the physical sheet and 
appears as a boundary bound state (see A. 2). Note that the derivation of the Cutkosky 
rules given above is perturbative, but as it rests on general principles such as analiticity and 
unitarity, we expect the rules to extend to singularities corresponding to such nonperturbative 
intermediate states as well, just as in bulk theories. In fact, the applications presented in the 
next Section require such an extension of Cutkosky rules to singularities with complex values 
of incoming momenta, but still in the physical sheet. 

Finally, let us remark that the derivation of the Landau equations and the Nakanishi proof 
for the Cutkosky rules for physical region singularities can be performed in a perturbation 
theory that starts around a nonzero value of A. In this case the reflected bulk propagator 
must be changed to 

i n r — i A 

kl — K 2 r — M 2 + ie K r + iX 

For amplitudes in the physical region, the last factor is just a nonsingular (phase-valued) 
expression, and so does not affect the derivations in any essential way. 

5 Applications 

In the following we restrict ourselves to 1 + 1 spacetime dimensions, since the example theory 
we use is sine-Gordon theory with integrable boundary conditions. The above formalism is 
of course more general, however, integrable field theories make possible an exact and non- 
perturbative verification of the general scheme. 
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5.1 Boundary bound state poles 



The simplest application of the Cutkosky rules is when there is a single intermediate line 
on-shell. If that is a bulk line, then we get back the well-known result for a first order pole 
corresponding to an intermediate particle in the bulk. If the line in question is a boundary 
line, we get a discontinuity of the form 

corresponding to the diagram Fig. 5.1, where E = Mcosh$ and m = Mcosu for a pole at 
imaginary rapidity ft = iu. 

This reproduces the result of Ghoshal and Zamolodchikov [20], which states that an inter- 
mediate on-shell boundary state produces a pole in the reflection factor R {&) of the form 



2-&-iu 



where the dimensionless couplings f\p are related to g\p by some normalization factors, coming 
partially from converting the 5^ to a function of d and also from relating the reflection factor 
in terms of the Feynman amplitude, but the precise form of this relation is not important for 
us here. 



5.2 Coleman-Thun diagrams 

We present some examples for more complicated intermediate diagrams following Coleman and 
Thun [1] (see also [2, 3, 22] for the boundary case). To be specific, we consider sine-Gordon 
theory with integrable boundary conditions. The groundwork for determining its spectrum 
and scattering amplitudes was laid down in [20, 21], while the bootstrap closure of this theory 
for general boundary conditions was obtained in [22], following the work in [3]. 

Before proceeding to the examples, we introduce a useful notation. Let /(#) be a mero- 
morphic function of one complex variable Let us suppose that the Laurent expansion of / 
around $ = i?o takes the form 

/(tf)=A(0-0 o ) n + X; M#-0o)* , A^O. 

k>n 

We then define 

/ [tfo] = A. 

If is not a pole or a zero of /, then this is just the value of the function. If $0 is a first 
order pole, then this is just the residue, while for a first order zero it is the derivative of / at 
00- 

The first family we consider are first order poles independent of the boundary conditions 
and occur in reflection factors of breathers on any boundary state (ground state and excited 
states alike). Let us denote the reflection factor of a breather B k with rapidity d by R^ (1?) 
(we omit the specification of the boundary state, as the argument does not depend on it). Then 
R^ has a pole at "& = inir/2X , n = 1, . . . ,k — 1. We pick a pole with a given n and write 
k = m + n. Following [22] this pole can be associated with the scattering process in diagram 
5.2 (we recall that an imaginary rapidity difference 1? = iu corresponds to a Euclidean angle u 
in the plane). The diagram obviously gives a first-order pole, since it has 3 propagators (one 
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Figure 5.2: Coleman-Thun graph for a boundary-independent pole in the breather reflection 
factor. The angles are uq = u\ = wk/2X and U2 = mix /2X. 




(a) u 7^ uq (b) u — > uq limit 



Figure 5.3: Calculating the contribution of the triangle graph using the bootstrap 
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Figure 5.4: Bootstrap relation for R( m+n ) 

5^ factors each), and a loop integral which means 2 integrations, leaving us with a single S 
function in the discontinuity which is characteristic of a first order pole. 

To establish that this process explains the pole of the reflection factor, we need to calculate 
its contribution to the residue. The straightforward method would be to write down the 
Cutkosky rules directly and calculate the integrals, which do give a contribution of the right 
form. However, the exact coefficient is not easy to obtain as one needs to ensure the right 
normalization of the reflection amplitude first, and then also extract the breather-breather 
fusion coupling <?^ m from the S matrix, once again taking care of all details of normalization. 

Instead, we can follow a simpler route, noting that the S matrices and the reflection 
amplitudes satisfy the bootstrap which entails a vast number of identities. We can tune the 
incoming rapidity away from $ = iuq, obtaining the process in Fig. 5.3 (a), which, by virtue 
of factorized scattering, has the amplitude 

gZ+ m R {n) (ti + i Ul )R {m) ($ - iu 2 ) S (2$ + i ( Ul - u 2 )) , ■& = iu 

Taking d — > iu 2 we hit a pole in the S matrix. The contribution of this pole reads 

g n n TR [n) (Ku 2 + m)) (-1)±S [i ( Ul + u 2 )\ 

On the other hand, in the same limit the diagram in Fig. 5.3 (a) becomes the one drawn in 
Fig. 5.3 (b), from which amplitude of the triangle graph Fig. 5.2 can be obtained by dividing 
with the breather fusion coupling gnt™ '• This means that to explain the pole in the reflection 
factor of B m+n , the following identity must hold: 

i?( m+n ) [iu 2 ] = (i(u 2 + Ul )) S [i (m + u 2 )\ 

which can be verified by direct substitution. 

It is very interesting to note that this identity is guaranteed by the bootstrap construction 
of the reflection factor R( m + n ) ; using the relation drawn in Fig. 5.4. 

We have also considered two other classes of Coleman-Thun diagrams, using similar argu- 
ments. Fig. 5.5 (a) shows a process where an incoming breather splits into two breathers of 
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B 



lk> 

(w_ r V/2 
ll> 



(a) Breather direct 



(b) Breather decay 



Figure 5.5: Coleman-Thun explanation for boundary dependent poles in breather reflection 
factors 



lower index, one of which excites the boundary state on which the other gets reflected. The 
process in Fig. 5.5 (b) describes an incoming boundary state decaying by the emission of a 
breather that fuses with the incoming breather to a higher one, which in turn gets reflected 
from the boundary. The notations for the boundary states \n) and the fusing angles v^, wi are 
explained in [18]. 

Naively, these diagrams would give a second order pole since the number of propagators is 
6 (including propagators for boundary states as dictated by the Cutkosky rules), while there 
are two loop integrals which give the order 6 — 2x2 = 2. However, the reflection at the 
middle of these diagrams takes place at a rapidity at which the corresponding amplitude has 
a first-order zero, which results in a first-order contribution. 

Straightforward application of bootstrap principles shows that the contribution of these 
diagrams to the given residues can be written as 



-R 



(n - l - k) [{w- k - *„_,) /2] A<£*> [(«,, - v k ) /2] S n - l ~ k ' l+k 



\k) 



/2A] 



(a): 

If these diagrams are to explain the residues, then the above expressions must be equal to 



XY K w -l - "*) /2] ^fft +fc "° [(^-« " w ~i) / 2 ] Sn,h ~ l [(« + *" 0vr/2A] 



(a): -!/„_,) /2] 



(6): 



We checked these equalities by direct evaluation. As before, they are also a consequence of 
the bootstrap relations between the reflection factors and S matrices. 

It is quite remarkable that the identities necessary for the diagrams to explain the residues 
of the pole hold by the bootstrap. This seems to be the case for breather reflection factors in 
sine-Gordon theory for generic values of the boundary parameters; although we could not check 
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it for all the poles in the reflection factors of all the breathers on generic excited boundary 
states due to the very tedious details, it was nevertheless true in each case we examined. For 
solitons, one does not expect this to be the case as in some diagrams the order of the pole is 
reduced by a cancellation between pairs of diagrams that are charge conjugate (in terms of 
soliton charge) analogous to the bulk Coleman-Thun diagrams, and so only the sub-leading 
term survives, which cannot be calculated so simply as above (in fact, such sub-leading terms 
have never been calculated in the literature, at least to our knowledge). 

There are also known cases when a certain pole in the reflection factor needs a combined 
explanation in terms of a sum of a Coleman-Thun diagram and a boundary bound state at the 
same time. For an example in the scaling Lee- Yang model cf. [2]. Such a situation also occurs 
in boundary sine-Gordon theory with Neumann boundary conditions, where it results from 
the confluence of two poles and a zero in the appropriate reflection factor as the boundary 
parameters are tuned to their value at the Neumann boundary condition, and has already 
been treated in detail in [23]. 

6 Conclusions 

The perturbative approach to boundary QFT around the Neumann boundary condition has 
some distinct advantages. It is very well suited for investigation of principles of boundary 
QFT, and for the extension of bulk results like LSZ reduction formulae, Landau equations 
and Cutkosky rules. For many applications it may be more convenient to start from a dif- 
ferent boundary condition (as is usually done in the literature), because that may reduce the 
number of diagrams (representing a partial re-summation of the perturbative series around 
the Neumann boundary condition) or eliminate tadpoles in case of a nontrivial vacuum field 
configuration (see Appendix A). 

The main result of this paper is the derivation of the Coleman-Norton interpretation, 
Cutkosky rules and their application to the boundary bootstrap. Although the rules are 
derived from perturbation theory, they are actually dependent only on some analyticity prop- 
erties and unitarity (this is especially obvious for Cutkosky rules, where the proof explicitly 
relies on unitarity). Therefore one expects that there is a non-perturbative formulation along 
the lines followed by Eden et al. in the last chapter of their book [4], which would be interest- 
ing to work out in detail. Such a formulation, besides being non-perturbative, would also give 
us an extension of the Cutkosky rules to the non-physical domain and to topologies distinct 
from the one considered in the proof of Section 4. 

Another open issue is to formulate a complete theory of perturbative renormalization with 
boundaries. In the perturbation theory around the Neumann boundary condition, it seems 
relatively easy to perform arguments similar to those in bulk theories, e.g. power counting 
analysis of divergent graphs. We only considered one-loop renormalization in boundary sine- 
Gordon theory in Appendix B, which was enough for a comparison with semiclassical results 
in the literature. It would be interesting to have a thorough investigation of renormalizability, 
proof of locality of counter-terms, investigation of possible anomalies of symmetries and so 
on. Such an investigation could provide e.g. a proof that the tadpoles indeed restore the 
correct vacuum expectation value of the field in the generic case (as shown for the toy model 
in Appendix A) thereby strengthening the derivation of analytic results from perturbation 
theory. 

Finally let us remark that Dirichlet boundary conditions cannot be reached by perturba- 
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tion theory from the Neumann this would require taking some coupling constants to 

infinity. However, it is rather trivial to extend the formalism of the paper to this case, which 
only requires the inclusion of a minus sign in the reflected propagator for fields satisfying 
Dirichlet boundary conditions. One also needs to require that the corresponding asymptotic 
fields satisfy Dirichlet boundary conditions, since if one insisted on an asymptotic Neumann 
boundary condition, that would mean adiabatic switching off for an infinitely strong inter- 
action, which would be a rather odd notion. Instead, taking the asymptotic field to satisfy 
Dirichlet boundary condition enables one to carry all the formalism over to this case. We 
only omitted this from the main text to maintain the arguments and the notation as simple 
as possible. 
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A A toy model 

In our approach, the free field has Neumann boundary conditions, and all the boundary in- 
teractions (including their quadratic part) are taken into account perturbatively. To illustrate 
how this procedure reproduces the interacting boundary conditions, we consider a toy model 
of a real scalar field <&. For simplicity, we take 1 + 1 spacetime dimensions (generalization 
to D + 1 is straightforward, see at the end of the section). We take the Lagrangian to be 
quadratic in the field 

r0 

dx 



/ 

J — oo 



\ ( w 2 - \ (W - ^ 



\®{y = o, tf 



so the model is exactly solvable. We remark that this toy model has already been used in the 
literature to illustrate several issues of boundary field theory [24]. 

A.l Exact solution 

The equations of motion are 

-d^ + d 2 y ^-m 2 ^ = , d y <f>\ y=0 = -\<f>{y = 0,t) (A.l) 
with the following solution for the modes of the field: 

fjy) = 1 ((k + i\)e iKy + (k - i\)e- iKy ) . 

For A < 0, there also exists a normalizable mode 

f B (y) = V2\\\e- X y 

The field has the mode expansion 



*(y, t) = Mt)f B {y)+ / 

J o 



dK4> K (t)f K {y) 
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where dynamics of the modes is 



dt 2 

d 2 Mt) 
dt 2 



= -(re 2 +m 2 )<M*) 
= -(m 2 -\ 2 )Mt) 



We remark that periodicity of the boundary mode in time requires A 2 < m 2 i.e. —m < A < 
(otherwise the vacuum $ = is unstable under perturbations in the direction of the boundary 
mode). The Lagrangian can be written in terms of modes as 



L = \ /°° dK (^ (t)2 " (fc2 + m2 )^W 2 ) + \ (^ (t)2 " n2 Mt) 2 ) , 

We can quantize the theory by introducing the conjugate momenta 

7r K (i) = 4> K {t) , n B (t) = (}>B(t) 
Creation/annihilation operators can be introduced: 

2lo^ \ J 



A 2 



o B {l) - -4= ( be~ int + &V n *) 



and satisfy 



a(re), <r («/) = <5(re — re') , 
For the case A > (i.e. no bound states) the propagator takes the form 
<0| T (d>(y, t)<f>(y', t>)) |0> = 9(t - t') <0| <f(y, *)<&(*/', |0) + 0(i' - t) <0| <&(*/', i) |0) 



oo />oo 



dudK 



oo — oo 



(27r) 2 uj 2 — k 2 — m 2 + ie 



g iK(y-y')-iu>(t-t') 



+ 



re — zA 



,iK-{y+y')-iuj(t-t') 



K + i\ I 

from which we can read off the reflection factor 

^/ n k — i\ 

W K + lX 

which has a bound state pole at 

k = msinhi? = — iX 
where d is the rapidity. Introducing u = —I'd 



A 



sin u = — - 



m 



For a bound state pole in the physical strip < u < ^ we need — m < A < 0. When this is 
satisfied, we get an additional term in the propagator. This can be obtained by continuing A 
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to negative values or by explicitely adding the contribution of the bound state mode, both of 
which gives 

{0\T^(y,tMy',t'))\0) = 



+ 

The bound state term can also be written in the form 

which means that (apart from the y dependence) the boundary bound state behaves as a free 
boundary field of mass Q propagating in + 1 dimensions (i.e. a harmonic oscillator) just as 
we postulated for the boundary degrees of freedom in our general exposition of perturbation 
theory, except that here the boundary degree of freedom arises as a boundary bound state and 
not from a separate field introduced in the Lagrangian. As a result, it is not sharply localized 
to the boundary, but its contribution decreases exponentially away from y = 0. 

A. 2 Boundary perturbation theory 

The reflection factor can be expanded in the coupling A 

R(k) = 1 + 2'£(- % -±\ (A.3) 

71=1 ^ ' 

We now proceed to show that this is correctly obtained using the perturbation theory intro- 
duced above. We suppose that A > 0, so that no boundary bound state exists, which would 
be a non-perturbative effect which perturbation theory is not expected to reproduce. 

A. 2.1 First order correction 

The interacting propagator can be written as 

G(y,t;y',t') = <0| T (<&(*/, t)<f>(y>, t')) exp j-^ dr : $(0, rf X |0> 

00 i / i\\ n n r°° 

the first order correction is 

i A f°° 

GW(y, t; y', t') = -- J dr (0| T ($(y, 0) = $(0,r) 2 : |0> 

/oo 
drG(y, t, 0, r)G(y', t', 0, r) 
-oo 



9(t - f) (0| $(y, my', t') |0) + 9(t' - t) (0| $(y, t)<f>(y', t') |0> 

f°° f°° dwdK i ( iK (y-y>)-Mt-n 

+ K^iX e i K (y + y')- iu} (t-t')\ 

k + iX J 

\^ e -\(y+y') (t _ t ') e -m-t) + 9(t>- t) e -^(*'-*)) (A.2) 
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Figure A.l: First order Feynman graphs for the propagator of the toy model 

where 

is the sum of the direct and reflected propagators. This contribution corresponds to the 
Feynman graphs in Fig. A.l. We can integrate r out to obtain 

r dT °^ t-WV' = 4 / ^ ^ 2 2 - 2 T ■ 2 / 2 T ■ e'^'-Mt-t') 

,/_ 00 7 (27r) d cj z — k z — m 2 + it uj z — k' z — m z + le 

The integral over k' can be performed using the residue theorem; since y' < 0, the contour 
must be closed in the lower half plane Siiik' < 0. Some further manipulation yields 



f dnduj 1 i 
J (2^^w 2 -k 2 - 



which means that the reflection factor becomes, at this order 



e iK(y+y')-iu(t-t') 



= 1 

K 



which agrees with (A. 3). 



A. 2. 2 Summing up to all orders 

At the nth order we get the contribution 



G^(y,t;y',t') = -^j (0| T ($(y, t)$(y', t')) JJ J dr t : <S>(0,n) 2 : |0) 

The relevant graphs are of the form shown in Fig. A. 2 . The number of contractions can be 
easily calculated: first one has to decide the order in which the vertices are to be connected 
starting from one end to the other: this can be done n! ways. Since each vertex has two 
identical legs, one has an additional factor of 2 n . Collecting all the terms we have 

G^(y,t;y',t') = ^ ("yj J Jn ■ • • J Jr n G(y, t - n ) 
x (jlG^n+i-Ti^Giy^Tn-t') 

/oo poo ( n ~y \ 

dn . . . / dr n G(y, t - n) H G (°> T ^ ~ n) G(y f , r n - t') 
-oo J— oo \j=l / 
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' 1 

y,t ' 

Figure A. 2: A typical diagram of order X n . To obtain all the diagrams, every line should be 
allowed to be either a direct bulk (solid) or a reflected bulk (dashed) line, which results in 
2 n+1 terms. 

The propagators G are the sum of the direct and reflected piece, so the formula contains 2 n+1 
terms, which turn out to give identical contributions in this case. The r integrals can be 
performed one by one, using the same method as for the first-order contribution, with the 
result: 

G (B) (V, t; y>, t>)= 2^ (- l ±Y 2 / e *WW) (A.4) 

J (2tt) z \ 2k J lo z — k z — m z + ie 

Summing up all the contributions, we obtain 

oo 

G(y,t;y',t') = £gW(j/, t;y',t') 

n=0 

f°° f°° dujdK * ( iK(y-y')-iu,(t-t') 

L^Loo (27r)2 W 2_ K 2_ m 2 + ie ^ 

+ R(k) e ^(y+y')-Mt-t')^ 

where 

oo / -\\n \ 
n=l 

which is exactly the result we expected. This means that the interacting boundary field, built 
by summing up the perturbative expansion around the free field with Neumann boundary 
conditions, does indeed satisfy the correct boundary conditions (A.l). 

We remark that this result can be continued analytically to A < 0: the pole at k = iX 
then crosses the contour of the k integral, and we obtain the result (A.l), which includes the 
contribution from the boundary bound state. So in this model, although the boundary bound 
state cannot be obtained by perturbation theory, it can be obtained by analytically continuing 
the resummed perturbative result from the regime where there is boundary bound state. It 
is important to resum the perturbation series: to any finite order the nontrivial singularity at 
k = —iX in the reflection factor is absent. 

A. 2. 3 Concluding remarks 

From (A.4) it is apparent that there is an infrared divergence at k = 0. This is a general 
feature of the perturbation theory around the Neumann boundary condition. The perturbative 
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expansion is a series in X/k and so it is convergent only for k > A. For long wavelength modes 
the expansion must be summed up and analytically continued. The physical manifestation of 
this problem is that R(k = 0) = —1, while the Neumann reflection factor is identically +1 for 
any k. The effective strength of the boundary interaction increases with the wavelength. 

To calculate the coordinate space two-point function (A. 4) one needs to integrate over all 
wavelengths and this can only be performed after summing up to all orders. In fact, the proper 
expression for the two-point function must include an infrared regulator, which can be removed 
only after summing up the leading behaviour of the function around n = 0. As an alternative, 
one can perform the whole calculation in momentum space, where every contribution is finite. 

For any finite order in perturbation theory, the amplitudes display a singularity at n = 0, 
which is in fact a solution to the Landau equations of Section 3.1 with all the bulk internal 
lines being on-mass-shell, and can be accounted for using the general formalism. Summing up 
the perturbation series moves this singularity to n = —iX, and in the regime — m < A < it 
describes a boundary excited state. 

Finally we note that the whole analysis can be generalized to D + 1 dimensions. The 
D + 1-momentum can be written as (u, k, n) and calculations for any mode with a given value 
of k are isomorphic to a calculation in the 1 + 1 dimensional case, with the particle mass m 
replaced by + m 2 . 



A. 3 Extension of the toy model: background fields 

We can extend the toy model to include a linear coupling at the boundary. This means that 
there is a classical background field in the vacuum, which poses an interesting question, since 
our perturbation theory works in an expansion around the Neumann boundary condition, 
which has no such fields. The question is whether with an appropriate re-summation in 
the boundary coupling we can get back to the correct value of the background field. Our 
Lagrangian is 



J — c 



f0 

L = I dx 

-oo 



-a$(j/ = 0,t)-^(y = 0, tf 



while the equations of motion take the form 

+ SjS - m 2 $ = , d y $\ y=0 = -a - A$(y = 0, *) (A.5) 

which has the classical vacuum solution 

cj> vac (x,t) = ^—e mx (A.6) 

m + A 

Let us try to compute this in perturbation theory in both a and A, i.e. around the Neumann 
boundary condition. The diagrammatic expression is 
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where the diagrams must be understood as classes in which every bulk line can be a direct 
or a reflected line (see Section 4). There is only a first order contribution in a, while the sum 
of all orders in A can be easily seen to give the result of the previous section for the two-point 
function, with a single a vertex attached to it and integrated over the boundary: 

/+oo /-+0O J r+oo J ■ 
dt — - = — = e~ iujt \e iKX + RU)e~ iKX \ 

R(k) = ^ 

The t and u> integrals can be trivially performed, and amount to substituting u = in the 
integrand, while the k integral can be easily performed using the residue theorem with the 
final result 

<<&(*,*)> " ""■ 



m + A 

which is fully consistent with (A. 6). 

Obviously, in a general theory if there is a nontrivial vacuum solution to the equations of 
motion, this means that we are going to have nontrivial tadpole diagrams in the perturbation 
expansion around the free field with Neumann boundary condition. It would be possible to 
eliminate tadpoles by expanding in fluctuations around this classical background field, but 
that would mean space dependent bulk couplings, which would make analytic investigation of 
perturbation theory extremely difficult. Therefore we choose to expand around the Neumann 
boundary condition, since re-summation of tadpoles should give us the same result as the 
background field method, as illustrated in the above example. In the toy model it can be easily 
seen (just by drawing all diagrams and using the above result for the one-point function) that 
all the higher correlation functions also get the correct tadpole contributions from perturbation 
theory. 



B One loop renormalization in sine-Gordon theory 

In [17] Kormos and Palla considered the semiclassical quantization of the two lowest energy 
static solutions of the boundary sine-Gordon model: 

m 2 r° 3 
V = Jij (l-cos/0$) ; *7 = M (l-cos|(<I>-v?o)) 

In their work the semiclassical energy corrections are obtained by summing up the contribu- 
tions of the oscillators associated to the linearized fluctuations around the static solutions. 
The appearing standard UV divergences, which are due to the non normal ordered nature 
of the Lagrangian, are canceled by counter terms: the coupling constants are renormalized 
m 2 — ► m 2 + 5m 2 and Mo — ► Mq + 5Mq. 5m 2 is chosen the same as in the bulk theory (obtained 
in standard perturbation theory), since due to the local nature of the counter term it cannot 
depend on the presence of the boundary. To determine 5Mq they impose the cancellation of 
logarithmic divergences. As a result the renormalized quantities of order 1 take the following 
form 



4vr J Q yjk 2 + m 2 ' 87r J yjk 2 + 



m 



where A is the momentum cutoff. 
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<k,(0) 



(a) 



(b) 



Figure B.l: Bulk divergence (a) and counter-term (b) at one loop 



Similar considerations lead to the same counter-terms in the analytically continued sinh- 
Gordon theory [15, 16]. The aim of this appendix is to derive these conjectured formulae in 
our systematic perturbative framework. We expand the Lagrangian to the appropriate order 
in/3 2 : 



V 



m 
T 



m 2 (3 2 



■$ 4 



U 



<I> /J 



4! " ' " 2 • 2 2 " 4! • 2 4 

where for simplicity we take cpo = 0. The divergent term at order m 2 (3 2 comes from diagram 
B.l (a), which is regularized by a momentum cutoff, and is canceled by a counter-term of the 
form 

ir 5m 2 2 
V C T = — ® 

(diagram B.l (b)). 

To compute 5m 2 we observe that on diagram B.l (a) the momentum (k,u) is incoming and 
outgoing in the same time, that is it does not contribute to the delta function: the integration 
is unconstrained (as opposed to the case where the loop propagator is changed to the reflected 
one, which results in a diagram that is not divergent at all). As a result the delta function 
factor is the same for both diagrams and from the cancellation of the vertex contribution we 
have 

5m 2 , „ -m 2 (3 2 i K dk /*°° duj i 
4! J -a 2^ J-c 



-4-3 



2tt uj 2 — k 2 — m 2 + ie 



(with the combinatorial factors given explicitly). Performing the uj integration we recover 
the standard result. We have three other divergent diagrams in this order, which can be 
obtained from B.l (a) by changing any or both of its external legs to the reflected propagator. 
Performing the same changes on diagram B.l (b), however, results a diagram which removes 
the required divergence. 

The divergent diagrams originating from the boundary term are of order Mo/3 4 and are 
canceled by a counter term of the form 



UcT 



5M p 2 
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(a) (b) (c) 



Figure B.2: Divergent contributions (a,b) and counter term (c) for the boundary interaction 

The related Feynman graphs are presented on diagram B.2 (a), (b) and (c). 

There are two differences compared to the bulk vertices. First of all, (5 has been replaced 
by ^ . More importantly, there is no momentum conservation in the boundary vertices and 
as a consequence not only diagram (a) but also diagram (b) is divergent, moreover they have 
the same contribution. Summing up these two terms the counter-term acquires a factor 2 (in 
addition to (3 — > ^) compared to the bulk computation and confirms the result (B.l). 
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